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Abstract Charm mesons produced in nucleus-nucleus collisions are expected to be less attenuated 
(quenched) by the medium than hadrons containing only light quarks, since radiative energy loss of heavy 
quarks should be reduced by the 'dead-cone' effect. We start from a published energy-loss model to derive 
the quenching for D mesons at the LHC, introducing an approximation of the dead-cone effect and employ- 
ing a Glauber-based description of the geometry of central Pb-Pb collisions to estimate the in-medium 
path lengths of c quarks. We show that the exclusive reconstruction of D° K~tt^ decays in ALICE 
allows to measure the nuclear modification factor of the D mesons transverse momentum distribution and 
the D/ charged hadrons ratio and, thus, to investigate the energy loss of c quarks. 



1 Introduction 



The ALICE experiment PP at the LHC will study nucleus- 
nucleus (AA) collisions at a centre-of-mass energy ^Snn = 
5.5 TeV (for Pb-Pb) per nucleon-nucleon (NN) pair in or- 
der to investigate the properties of QCD matter at energy 
densities of up to several hundred times the density of 
atomic nuclei. In these conditions a dcconfined state of 
quarks and gluons is expected to be formed. 

Hard partons and heavy quarks, abundantly produced 
at LHC energies in initial hard-scattering processes, are 
sensitive probes of the medium formed in the collision 
as they may lose energy by gluon bremsstrahlung while 
propagating through the medium j^lEMlE]- The attenua- 
tion (quenching) of leading hadrons and jets observed at 
RHIC IHHZIIHlini is thought to be due to such a mecha- 
nism. The large masses of the charm and beauty quarks 
make them qualitatively different probes, since, on well- 
established QCD grounds, in-medium energy loss of mas- 
sive partons is expected to be significantly smaller than 
that of 'massless' partons (light quarks and gluons) jlO) . 
Therefore, a comparative study of the attenuation of mass- 
less and massive probes is a promising tool to test the 
consistency of the interpretation of quenching effects as 
energy loss in a deconfined medium and to further inves- 
tigate the properties (density) of such medium. 

In the first part of this paper, we shortly summarize 
a widely used model of parton energy loss and we discuss 
how we apply it in our simulation. In the second part, 
we show that the exclusive reconstruction of D° K~7t~^ 
decays with ALICE allows to carry out the mentioned 
comparative quenching studies by measuring: 



the nuclear modification factor of D mesons as a func- 
tion of transverse momentum (pt) 



_ dA^AA/dpt /binary NN collision 

i?AA(Pt = , 

dNpp/dpt 



(1) 



which would be equal to 1 if the AA collision was a 
mere superposition of independent NN collisions, with- 
out nuclear or medium effects; in central Au-Au col- 
lisions at RHIC (^i^ 200 GeV) Raa ^ 0.2 for 
both 7r° and charged hadrons in the range 4 < pt < 
8 GeV/c EHmHI; 

the ratio of the nuclear modification factors of D mesons 
and of charged (non-charm) hadrons, as a function of 
Pf- 

(2) 



i?D/^(pt) = i?2Abt)/i?^A(Pt); 



hereafter, this quantity is called D/ charged hadrons 
(D/h) ratio. 



2 Parton energy loss and the dead-cone 
effect for heavy quarks 

In this work we use the quenching probabilities (or weights) 
calculated in Ref. in the framework of the 'BDMPS' 
(Baier-Dokshitzer-Mueller-Peigne-Schiff) formalism sum- 
marized in the following. The energy loss obtained with 
the quenching weights is presented in Section |21 

An energetic parton produced in a hard collision ra- 
diates a gluon with a probability proportional to its path 
length L in the dense medium. Then, the radiated gluon 
undergoes multiple scatterings in the medium, in a Brown- 
ian-like motion with mean free path A which decreases as 
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the density of the medium increases. The number of scat- 
terings of the radiated gluon is also proportional to L. 
Therefore, the average energy loss of the parton is pro- 
portional to i^. 

The scale of the energy loss is set by the 'maximum' 
energy of the radiated gluons, which depends on L and on 
the properties of the medium: 



(3) 



where q is the transport coefficient of the medium, defined 
as the average transverse momentum squared transferred 
to the projectile per unit path length, q = (^t )inedium / ^ C3 

In the case of a static medium, the distribution of the 
energy w of the radiated gluons (for u <C liJc) is of the 
form: 

_ d/ ^ 2 as Cr fuj^ 

did TT Y 2w ' 

where Cr is the QCD coupling factor (Casimir factor), 
equal to 4/3 for quark-gluon coupling and to 3 for gluon- 
gluon coupling. The integral of the energy distribution up 
to Wc estimates the average energy loss of the parton: 



{AE) 



dl 



Lo—duo ca agCRUc cti oisCRqlF' . (5) 
dcj 



The average energy loss is: proportional to as Cr and, 
thus, larger by a factor 9/4 = 2.25 for gluons than for 
quarks; proportional to the transport coefficient of the 
medium; proportional to L^; independent of the parton 
initial energy E. It is a general feature of all parton en- 
ergy loss calculations P|r^|IHt ^lllll2j that the energy dis- 
tribution ^ docs not depend on E. Depending on how 
the kinematic bounds are taken into account, the result- 
ing AE is £'-independent [H| or depends logarithmically 
on E |12| . However, there is always a stronger intrin- 
sic dependence of the radiated energy on the initial en- 
ergy, determined by the fact that the former cannot be 
larger than the latter, AE < E. Within the above toy- 
model derivation which agrees with the main features of 
the BDMPS formalism, this kinematic constraint could 
be included by truncating the gluon energy distribution 
uidl/doj at min(a;c,£') rather than at its natural upper 
limit iOc- This would give, from Q and ©, {AE) oc 
as Cr ^/oj^ ■\/™i^('-^cj E). For E < o^c, we have {AE) oc 
y/qy/E L: the kinematic constraint turns the L-dependence 
from quadratic to linear. As we shall discuss in Section |S1 
the kinematic constraint can be equivalently interpreted 
as a reduction of the effective path length in the medium, 
i.e. of the length along which the energy of the parton is 
larger than zero and gluons can be radiated. 

The transport coefficient is proportional to the den- 
sity of the scattering centres and to the typical momen- 
tum transfer in the gluon scattering off these centres. A 
review of the estimates for the value of the transport co- 
efficient in media of different densities can be found in 
Ref. : the estimate is gcoid — 0.05 GeV^/fm for cold 
nuclear matter and up to 0.5 GeV^/fm for a hadron gas; 
for a QGP formed at the LHC with energy density e ~ 50- 
100 GeV/fm^ q is expected to be of ~ 5-10 GcVVfm. 



The medium-induced energy loss of heavy quarks was 
first studied in Refs. |14lll5j . Later, in Ref. ^01 it was 
argued that for heavy quarks, because of their large mass, 
the radiative energy loss should be lower than for light 
quarks. The predicted consequence of this effect was an 
enhancement of the ratio of D mesons to pions (or hadrons 
in general) at moderately-large (5-10 GeV/c) transverse 
momenta, with respect to that observed in the absence of 
energy loss. 

Heavy quarks with momenta up to 40-50 GcV/c prop- 
agate with a velocity which is significantly smaller than 
the velocity of light. As a consequence, in the vacuum, 
gluon radiation at angles smaller than the ratio of their 
mass to their energy 0q = m/E is suppressed by destruc- 
tive quantum interference |16j . The relatively depopulated 
cone around the heavy-quark direction with 6* < 0o is 
called 'dead cone'. 

In Ref. the dead-cone effect is assumed to char- 
acterize also in-medium gluon radiation and the energy 
distribution of the radiated gluons for heavy quarks, 
is estimated to be suppressed by the factor: 



d/ 

doj 



Heavy/ 



'd/ 

duj 



Light 




■-1 -2 



(6) 



Fn/i.{m,E,uj,q) 



where the expression for the characteristic gluon emission 
angle O ~ {q / to'^)'^ ^ ^ has been used. The hcavy-to- 
light suppression factor i^H/L in I© increases (less sup- 
pression) as the heavy-quark energy E increases (the mass 
becomes negligible) and it decreases at large w, indicating 
that the high-energy part of the gluon radiation spectrum 
is drastically suppressed by the dead-cone effect. 

The first (indirect) measurements of charm produc- 
tion, in the semi-electronic decay channel, in Au-Au col- 
lisions at RHIC ^7] indicate (within the still large ex- 
perimental errors) that there may be no attenuation of 
the D meson pt distribution, contrarily to what observed 
for pions and charged hadrons. This result supports the 
scenario proposed in Ref. and it has stimulated con- 
siderable theoretical interest in the subject |18II19| . Our 
approach, described in the next section, is to implement 
in the energy-loss simulation an algorithm to account for 
the dead-cone effect. 



3 Simulation of energy loss 

The quenching weight |2()II1H is defined as the probability 
that a hard parton radiates an energy AE due to scat- 
tering in spatially-extended QCD matter. In Ref. ^l], the 
weights are calculated on the basis of the BDMPS formal- 
ism, taking into account both the finite in-medium path 
length L and the dynamic expansion of the medium after 
the nucleus-nucleus collision. The input parameters for 
the calculation are the length L, the transport coefficient 
q and the parton species (light quark or gluon). 
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Figure 1. Distribution of the path lengths in the transverse 
plane for partons produced in Pb-Pb collisions with impact 
parameter b < 3.5 fm. 



The distribution of the in-medium path length in the 
plane transverse to the beam line^ for central Pb-Pb col- 
lisions (impact parameter b < 3.5 fm. corresponding to 
the 5% most central collisions) is calculated in the frame- 
work of the Glauber model for the collision geometry j^J . 
For a given impact parameter, hard-parton production 
points are sampled according to the density Pcoii{x,y) of 
binary nucleon-nucleon collisions in the transverse plane 
and their azimuthal propagation directions are sampled 
uniformly. For a parton with production point (xq, yo) and 
azimuthal direction {ux, Uy), the path length is defined as: 



d; / pcoU (a;o + lux,Vi3 + l Uy) 
0.5 dlpco\i[xo + lux,yo + luy) 



(7) 



Many sampling iterations arc performed varying the 
impact parameter h from 0.25 fm to 3.25 fm in steps of 
0.5 fm. The obtained distributions are given a weight b, 
since we verified that dcr'^^'^'^ /d6 cx 5 for 6 < 3.5 fm, and 
summed. The result is shown in Fig.^ The average length 
is 4.5 fm, corresponding to about 70% of the radius of 
a Pb nucleus and the distribution is significantly shifted 
towards low values of L because a large fraction of the 
partons are produced in the periphery of the superposition 
region of the two nuclei ('corona' effect). 

The definition of L in lO is exact in the case of a 
cylindrical collision density profile (i.e. pcoii = Po 0{R — 
+ y"^) for a cylinder with radius R). It is approxi- 
mated for a realistic profile (derived from a Wood-Saxon 
nuclear density profile). The inclusion of nuclear geometry 
effects is not straight-forward in the scheme of the quench- 
ing weights, where the transport coefficient and the path 
length are considered as two distinct input parameters. 



^ Partons produced at central rapidities propagate in the 
transverse plane. 



Indeed, the lower medium density in the 'corona' could be 
modeled with a reduction of the transport coefficient or of 
the path length. However, since the energy loss is deter- 
mined by the quantities tOc Q L'^ and lOc L oa q , 
the optimal solution would be to include the geometry 
profile directly in the calculation of these two parameters, 
defining them as line integrals in dl similar to the integrals 
that appear in Q . For this work we adopt Q as a practi- 
cal definition of length in the medium, which allows us to 
employ a distribution of L rather than a constant value; in 
the following we show that this point is quite relevant. A 
more refined treatment of nuclear geometry is discussed in 
Ref. 1221, where, however, energy loss is modeled as a sim- 
ple exponential absorption that cannot be directly related 
to the medium properties. 

For a given value of the transport coefficient q and a 
given parton species, we use the numerical routine pro- 
vided in Ref. ^I] to calculate the energy-loss probability 
distribution P{AE; L) for the integer values of L up to 
15 fm. Then, these 15 distributions are weighted accord- 
ing to the path-length probability in Fig. ^ and summed 
up to obtain a global energy-loss probability distribution 
P{AE). The energy loss to be used for the quenching sim- 
ulation can be directly sampled from the P{AE) distri- 
bution corresponding to the chosen q and to the correct 
parton species. Figure |21 (left) reports P{AE) for light 
quarks and for gluons, as obtained with different values 
of the transport coefficient; in the figure, the 'peak' at 
AE = represents the probability to have no medium- 
induced gluon radiation. 

The predicted lower energy loss for charm quarks is 
accounted for by multiplying the P{AE) distribution for 
light quarks by the dead-cone suppression factor i^H/L 
given in with lo = AE. Since i^H/L depends on the 
heavy-quark energy E, the multiplication has to be done 
for each c quark or, more conveniently, in bins of the quark 
energy. It was verified that this multiplication is equiva- 
lent to recalculating the quenching weights with the gluon 
energy distribution for heavy quarks modified according 
to The results on the energy-loss probability distri- 
bution obtained with the two methods, 'Fh/l ® P{AE) 
product' and 'i^n/L in the quenching weights', were com- 
pared for g = 1 GeV^/fm, L = 5 fm, mc = 1.2 GeV/c^ 
and pt = 10, 20, 30 GeV/c |2H|- The comparison for 
Pt = 10 GeV/c, shown in Fig. El (right), is quite satis- 
factory. A similar agreement is found for pt = 20 and 
30 GeV/c. 

Figure |21 shows the average relative energy loss as a 
function of the transverse momentum for gluons, light 
quarks and charm quarks {nic = 1.2 GcV/c^). Differ- 
ent values of the transport coefficient q are considered in 
the various panels and the described dead-cone correction 
(energy-dependent Fh/l ^ P{AE) product) is used for 
charm quarks. For each value of the parton pt many values 
of AE are sampled according to the energy-loss probabil- 
ity distribution and the kinematic constraint AE < E, 
discussed in Section 12 is applied as follows: if AE > pt, 
then AE = pf The ratio of the relative energy losses 
for gluons and for light quarks is compatible with the ra- 
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Figure 2. Left-liand panel: energy-loss probability distribution, P{AE), for different values of q, for light quarks and for gluons 
(thicker dashed line, only for 4 GeV'^/fm). Right-hand panel: P{AE) for c quarks (rric = 1.2 GeV/c^, pt = 10 GeV/c) without 
dead cone and with two different implementations of the dead-cone effect (see text). 



tio of their Casimir factors (2.25) at high pt, where the 
kinematic constraint is not relevant. The ratio of the rel- 
ative energy losses for light quarks and c quarks increases 
with increasing q, particularly in the high-pt region, show- 
ing that the dcad-conc effect is medium-dependent. With 
q = 4 GeVVfm, our estimated transport coefficient for 
the LHC (see next paragraph), the average relative en- 
ergy loss is ~ (85 — 0.6pt(GcV/c))% for gluons, ~ (75 — 
0.8pt(GeV/c))% for light quarks and « 25-30% for c 
quarks. 

Remarkably, for charm quarks wc find that, for given g, 
the average relative energy loss is approximately indepen- 
dent of the quark energy, {AE/E) sa const., (see Fig. 
while for massless partons the average relative energy loss 
is clearly decreasing as the parton energy increases (the 
BDMPS average relative energy loss for massless partons 
would be {AE/E) oc 1/E] this dependence is, then, weak- 
ened by the kinematic constraint). On the basis of this 
observation we expect that not only the magnitude but 
also the pt-dependence of the nuclear modification factor 
i?AA of D mesons can be significantly affected by the dead 
cone. 

Before describing the estimation of the transport coef- 
ficient to be used in the simulation, we shortly comment on 
the implementation of the kinematic constraint. Taking, 
as we do, AE = pt when the sampled AE is larger than pt 
is equivalent to truncating the energy-loss probability dis- 
tribution P{AE) at AE = Pt and adding the (5-function 
6{AE - Pt) d£ P{£) to it. The total integral of P is, 
in this way, maintained equal to 1. The dependence of the 
final result on the kinematic constraint was argued |23l 
[TT] to illustrate the theoretical uncertainties. These can 
be sizeable (more than a factor 2) for low pt and suffi- 
ciently large q and L in the present study. They should 



be improved in a refined calculation but we did not try to 
quantify them in detail^. 

For the estimation of the transport coefficient q we 
require for central nucleus-nucleus collisions at the LHC 
a leading-particle quenching at least of the same mag- 
nitude as that observed at RHIC. We, therefore, derive 
the nuclear modification factor i?AA for charged hadrons 
produced at the LHC and we choose the transport co- 
efficient in order to obtain i?AA — 0.2-0.3 in the range 
5 < Pt < 10 GeV/c (for RHIC results see e.g. Refs. 

■SI)- 

The transverse momentum distributions, forpt > 5 GeV/c, 
of charged hadrons are generated by means of the chain: 

1. generation of a parton. quark or gluon, with pt > 
5 GeV/c, using PYTHIA [S^j proton-proton with s/s = 
5.5 TeV and CTEQ 4L parton distribution functions pS] : 
with these parameters, the parton composition given 
by PYTHIA is 78% gluons and 22% quarks (average 
over Pt > 5 GeV/c); 

2. sampling of an energy loss AE according to P{AE) 
and calculation of the new transverse momentum of 
the parton, p[ = pt — AE (if AE > pt, AE = pt and 

v't = 0); 

3. (independent) fragmentation of the parton to a hadron 
using the leading-order Kniehl-Kramer-Potter (KKP) 
fragmentation functions ,27,. 

Quenched and unqucnchcd pt distributions are obtained 
including or excluding the second step of the chain. Fig- 

^ These topics are considered also in Ref. 1241 of which the 
author became aware while finalizing the present work. Very 
useful discussions with the authors of Ref. |24| are acknowl- 
edged. 
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Figure 3. Average relative energy loss as function of the transverse momentum for gluons, light (massless) quarks and charm 
quarks {nic — 1.2 GeV/c^). 



ure 01 shows Raa for hadrons, calculated as the ratio of 
the Pt distribution with quenching to the pt distribution 
without quenching. Different values of q are considered 
in the left-hand panel of the figure: a value as large as 
4 GeV^/fm is necessary to have Raa — 0.25-0.3 in 5 < 
Pt < 10 GeV/c. In the right-hand panel, for g = 4 GeV^/fm, 
we compare the results obtained considering all partons 
as gluons or all partons as quarks, in order to remark and 
quantify the larger quenching of gluons with respect to 
quarks. 

Since the transport coefficient determines the size of 
the energy-loss effect, we shortly discuss the choice of 
9 = 4 GeV^/fm. This value corresponds, according to 
the estimates reported in Ref. ^31, to an energy den- 
sity e ~ 40-50 GeV/fm^, which is about a factor 2 lower 
than the maximum energy density expected for central 
Pb-Pb collisions at the LHC. The value is, therefore, rea- 
sonable. 



Using the same quenching weights, the suppression ob- 
served at RHIC is reproduced in Ref. ^Jj with the much 
lower value q = 0.75 GeV^/fm. However, the constant 
length L = 6 fm is there used, rather than a realistic dis- 
tribution of lengths, thus obtaining a significantly stronger 
quenching. This simple approximation captures the main 
features of Raa for 5-10% central collisions, which depend 
on the combination q (oc ujc L) more than on q and L 
separately )11| . In addition, transverse momentum distri- 
butions are steeper at RHIC energy than at LHC energy 
and, consequently, to obtain the same Raa suppression in 
the two cases one needs a larger q (or (AE)) at the LHC 
than at RHIC. 

Because of the kinematic constraint AE < E, the use 
of a constant length of the order of the nuclear radius or 
even the use of the average length from a detailed dis- 
tribution can produce quite different results with respect 
to those obtained taking into account the complete dis- 
tribution. This is demonstrated in Fig. [3 for q ~ 0.5- 
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Figure 4. Nuclear modification factor for charged hadrons for different values of q (left) and comparison of quark and gluon 
quenching for g = 4 GeV^/fm (right). 
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Figure 5. Nuclear modification factor for charged hadrons for different values of q as obtained using the complete Glauber-based 
distribution (from Fig. or constant values of L. 
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1 GcVVfm, L = 6 fm ffivcs almost a factor 2 difference in 
i?AA at pt ^ 10 GeV/c and the complete L distribution 
is equivalent to a constant length which decreases as q in- 
creases, 5, 4.5, 4, 3.5 fm for q = 0.2, 0.5, 1, 4 GeVVfm. 
This behaviour is clearly due to an upper 'cut-off' of the 
length distribution: large lengths correspond to very high 
values of AE, but, since AE cannot be higher than the 
initial parton energy £', large lengths arc not 'fully ex- 
ploited'; this corresponds to a cut-off; e.g. for many par- 
tons of moderate energy a length of 8 fm is equivalent 
to a length of 4 fm, because after propagating for 4 fm 
they have lost all their initial energy. As a consequence, 
the length distribution corresponds to an average 'effec- 
tive' length lower than its arithmetic average. The cut-off 
moves towards lower lengths as q increases and, thus, the 
average effective length decreases. For g = 4 GeV^/fm, the 
average effective energy loss is approximately: {AEes) ~ 
{AE) X (Lcff)^ / {Lf = {AE) X (3.5/4.5)2 w 0.6 {AE). 

Another important observation revealed by Fig. El is 
the fact that the use of the complete L distribution reduces 
the increase of i?AA with pt (at RHIC, i?AA is found to be 
independent of pt in the range 4-10 GeV/c [7|). This hap- 
pens because higher-energy partons can exploit the large- 
L tail more than lower-energy partons and, consequently, 
for them the cut-off is shifted towards larger lengths. 

Charm quarks are generated using PYTHIA, tuned in 
order to reproduce the single-inclusive c (and c) pt distri- 
bution predicted by the pQCD program HVQMNR gHI 
with m-c = 1.2 GeV/c^ and factorization and renormaliza- 
tion scales /iF = ^J'R = 2 mt = 2 + pi (the details on 
this tuning can be found in Ref. [SU]). In HVQMNR we use 
the CTEQ 5M parton distribution functions includ- 
ing, for Pb-Pb, the nuclear shadowing effect by means 
of the EKS98 parameterization and the parton in- 
trinsic transverse momentum broadening as reported in 
Ref. [nn|. With these parameters, the cc production yields 
in pp collisions at ^/s = 14 TeV and in central (5%) 
Pb-Pb coUisions at ^s^n = 5.5 TeV are iVpp = 0.16 and 

^Pb-Pb ~ 115' respectively. The yield given for Pb-Pb al- 
ready includes a 65% reduction due to shadowing. 

Energy loss for charm quarks is simulated following a 
slightly different procedure with respect to that for light 
quarks and gluons. Since the total number of cc pairs per 
event has to be conserved, in the cases where the sampled 
AE is larger than pt , we assume the c quark to be thermal- 
ized in the medium and we give it a transverse momentum 
according to the distribution dA^/dmt oc mt exp(— mt/T), 
as suggested in Ref. [TS|. We use T = 300 MeV as the ther- 
malization temperature for c quarks. The other difference 
with respect to the case of masslcss partons is that we 
use the standard string model in PYTHIA for the c-quark 
fragmentation (more details can be found in Ref. j33|). 



4 Charm reconstruction with ALICE 

The transverse momentum distribution of charm mesons 
produced at central rapidity, |y| < 1, can be directly mea- 
sured from the exclusive reconstruction of D*^ — s- K^tt"*" 



decays (and charge conjugates) in the Inner Tracking Sys- 
tem (ITS), Time Projection Chamber (TPC) and Time Of 
Flight (TOF) of the ALICE barrel (|?7| < 0.9). A schema of 
the employed detectors is reported in Fig. El The displaced 
vertices of D° decays (cr = 124 ^m) can be identified in 
the ITS with silicon pixels, that provide a measurement of 
the track impact parameters to the collision vertex with a 
resolution better than 50 /im for pt > 1 GcV/c. The low 
value of the magnetic field (0.4 T) and the K/tt separation 
in the TOF detector allow to extend the measurement of 
the D*^ production cross section down to almost zero trans- 
verse momentum. The strategy for this analysis and the 
selection cuts to be applied were studied with a realis- 
tic and detailed simulation of the detector geometry and 
response, including the main background sources j33ll34j . 

The expected performance for central Pb-Pb (6 < 
3.5 fm) at ^s^n — 5.5 TeV and pp collisions at ^/s ~ 
14 TeV is summarized in Fig. [3 The accessible pt range is 
1-14 GeV/c for Pb-Pb and 0.5-14 GeV/c for pp, without 
assuming dedicated triggers; triggering on high-pt tracks 
or on more specific kinematic and topological features us- 
ing the High Level Trigger may allow to further extend 
these Pt ranges. The statistical error corresponding to 1 
month of Pb-Pb data-taking and 9 months of pp data- 
taking is better than 15-20% and the systematic error 
(acceptance and efficiency corrections, subtraction of the 
feed-down from B D° -|- X decays, cross-section nor- 
malization, centrality selection for Pb-Pb) is better than 
20%. More details are given in Ref. |33| . 



5 Results: R^a and R-o/h 

The nuclear modification factor |^ for D*^ mesons is re- 
ported in Fig.lHl Nuclear shadowing, parton intrinsic trans- 
verse-momentum broadening and energy loss are included. 
The dead-cone effect is not included in the left-hand panel 
and included in right-hand panel. Different values of the 
transport coefficient are used for illustration; we remind 
that the value expected on the basis of the pion quenching 
observed at RHIC is g = 4 GeV^/fm. The reported sta- 
tistical (bars) and systematic (shaded area) errors are ob- 
tained combining the previously-mentioned errors in Pb-Pb 
and in pp collisions and considering that the contribu- 
tions due to cross-section normalization, feed-down from 
beauty decays and, partially, acceptance/efficiency correc- 
tions will cancel out in the ratio. An uncertainty of about 
5% introduced in the extrapolation of the pp results from 
14 TeV to 5.5 TeV by pQCD is also accounted for (see 
Ref. ESI). 

The effect of shadowing, clearly visible for q = (no 
energy loss) as a suppression of i?AA, is limited to pt < 6- 
7 GeV/c (using EKS98 EH)- Above this region, only pos- 
sible parton energy loss is expected to affect the nuclear 
modification factor of D mesons. 

For g = 4 GeV^/fm and no dead cone, we find i?AA 
reduced, with respect to 1, by a factor about 3 and slightly 
increasing withpt, from 0.3 at 6 GeV/c to 0.4 at 14 GeV/c. 
Even for a transport coefficient lower by a factor 4, g = 
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2 4 6 8 10 12 14 2 4 6 8 10 12 14 

Pt [GeV/c] p, [GeV/c] 

Figure 7. Double-differential cross section per nucleon-nucleon collision for D° production as a function of pt , as it can be 
measured with 10^ central Pb-Pb events (left), corresponding to 1-month data-taking, and 10^ pp minimum-bias events (right), 
corresponding to 9-months data-taking. Statistical (inner bars) and pt-dependent systematic errors (outer bars) are shown. A 
normalization error of 11% for Pb-Pb and 5% for pp is not shown. 



1 GeVVfm, Raa is significantly reduced (0.5-0.6). When 
the dead-cone effect is taken into account, the Raa re- 
duction due to quenching is found to be lower by about 
a factor 1.5-2.5, depending on q and pt- For our reference 
transport coefficient, 4 GeV^/fm, Raa with dead cone is 
equal to 0.6 and essentially flat as a function of pt. 

We point out that the estimated systematic uncer- 
tainty of about 18% may prevent from discriminating be- 
tween a scenario with moderate quenching and negligible 



dead-cone effect (e.g. q = 1 GeV^/fm in the left-hand 
panel of Fig. ^ and a scenario with large quenching but 
also strong dead-cone effect (e.g. q = A GeV'/fm in the 
right-hand panel). 

The comparison of the quenching of c-quark-originated 
mesons and massless-parton-originated hadrons will be 
the best-suited tool to disentangle the relative importance 
of energy- loss and dead-cone effects. T\ieT> / charged hadrons 
ratio i?D //, , defined as in (O , is presented in Fig. IHl for 
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Figure 8. Nuclear modification factor for D" mesons with shadowing, intrinsic kt broadening and parton energy 
panel: without dead cone correction; right-hand panel: with dead cone correction. Errors corresponding to the 
are shown: bars = statistical, shaded area = systematic. 
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Figure 9. Ratio of the nuclear modification factors for D" mesons and for charged hadrons. Left-hand panel: without dead cone 
correction; right-hand panel: with dead cone correction. Errors corresponding to the curve for q = 0.05 GeVVfm are shown: 
bars — statistical, shaded area — systematic. 



the range 5 < pt < 14 GeV/c. We used R^a calculated 
as previously described and R^j^, without and with dead 
cone, as reported in Fig. |H1 Being essentially a double ra- 
tio Pb-Pb/Pb-Pb X pp/pp, many of the systematic un- 
certainties on i?D//t cancel out (ccntrality selection and, 
partially, acceptance/efficiency corrections and energy ex- 
trapolation by pQCD). The residual systematic error is 
estimated to be of about 10-11%. 

We find that, if the dead-cone correction for c quarks 
is not included, Ru/h is essentially 1 in the considered pt 
range, independently of the value of the transport coeffi- 
cient, i.e. of the magnitude of the energy-loss effect. When 
the dead cone is taken into account, i?D//i is enhanced of 
a factor strongly dependent on the transport coefficient 
of the medium: e.g. 2-2.5 for (7 = 4 GcV'^/fm and 1.5 for 



q = 1 GeV^/fm. The enhancement is decreasing with pt, 
as expected (the c-quark mass becomes negligible). 

The Ru/h ratio is, therefore, found to be enhanced, 
with respect to 1, only by the dead cone and, consequently, 
it appears as a very clean tool to investigate and quantify 
this effect. 

Since hadrons come mainly from gluons while D mesons 
come from (c) quarks, the D/h ratio should, in principle, 
be enhanced also in absence of dead-cone effect, as a con- 
sequence of the larger energy loss of gluons with respect 
to quarks. Such enhancement is essentially not observed 
in the obtained Ru/h because it is 'compensated' by the 
harder fragmentation of charm quarks with respect to light 
quarks and, particularly, gluons. With z the typical mo- 
mentum fraction taken by the hadron in the fragmenta- 



10 



Andrea Dainese: Perspectives for the study of charm in-medium quenching at the LHC with ALICE 



hadron _ z p^'^^^"'^ , and AE the average energy loss 

r ,1 J / partonx/ parton 

tor the parton, (p^ Y — Pt 



tion, 



hadron \/ 



hadron 
Pt 



AE, we have 
zAE, 



(8) 



meaning that the energy loss observed in the nuclear modi- 
fication factor is, indeed, z AE. We have, thus, to compare 

— >D AEq to ^gluon — vhadron -^-^gluon- ^Vith ^gluon — *hadron ~ 

0.4, Zc^D ~ 0.8 forp°^'' > 5 GeV/c and AE^ = Z\£;giuon/2.25 
(without dead cone), we obtain 



luon — !-hadron 



AE, 



gluon ■ 



(9) 



This simple estimate confirms that the quenching for D 
mesons is almost the same as for (non-charm) hadrons, if 
the dead-cone effect is not considered. 

The errors reported in Fig. El show that ALICE is ex- 
pected to have good capabilities for the study of Ru/h- in 
the range 5 < pt < 10 GeV/c the enhancement due to the 
dead cone is an effect of more than 3 cr for g > 1 GeV^/fm. 
The comparison of the values for the transport coefficient 
extracted from the nuclear modification factor of charged 
hadrons and, independently, from the D / charged hadrons 
ratio can provide an important test for the coherence of 
our understanding of the energy loss of hard probes prop- 
agating in the dense QCD medium formed in Pb-Pb col- 
lisions at the LHC. 
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